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I— I ' Abstract 

^ : 

There are not many kinetic models where it is possible to prove bifurcation phenomena 
' ' for any value of the Knudsen number. Here we consider a binary mixture over a line with 

' collisions and long range repulsive interaction between different species. It undergoes a 

segregation phase transition at sufficiently low temperature. The spatially homogeneous 
I Maxwellian equilibrium corresponding to the mixed phase, minimizing the free energy at 

high temperature, changes into a maximizer when the temperature goes below a critical 
I I value, while non homogeneous minimizers, corresponding to coexisting segregated phases, 

^ ' arise. We prove that they are dynamically stable with respect to the Vlasov-Boltzmann 

T-H ' evolution, while the homogeneous equilibrium becomes dynamically unstable. 
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o^ : 

'. 1 Introduction and Main Results. 

> . 

^ . The phenomenon of phase transition in a thermodynamic system is usually described by the 

^ I arising of multiple minimizers of the free energy. Namely, when the temperature is lowered 

■ - - ' below a certain critical value, the unique equilibrium minimizer becomes a local maximizer 

and new minimizers appear. This is interpreted as a loss of stability of the old minimizer 
and the birth of new stable states. Here the meaning of stability is merely related to the free 
energy minimizing properties of the state. Of course it would be desirable to have a more 
detailed, dynamical analysis of the stability properties involved in such a phenomenon of phase 
transition. This issue is, generally speaking, very difficult, because the natural dynamics is 
a many-component one, whose detailed understanding is very far from being attainable. A 
simpler but, in our opinion, still interesting problem is to study this kind of behavior at 
kinetic level, where, instead of a system of a huge number of particles, one has to study a 
partial differential equation for the probability distribution of particles on the one particle 
phase space. This is the problem we want to address in this paper. 

The standard kinetic model of a rarefied gas undergoing collisions (Boltzmann equation), 
however, describes essentially the ideal gas, which does not exhibit phase transitions. Fol- 
lowing van der Waals, it is convenient to include a long range attractive interaction between 
particles to see a vapor-liquid transition. With the introduction of such an interaction, new 
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problems arise, due to the fact that nothing prevents the system from collapsing (Statistical 
Mechanics instabihty). The van der Waals approach of adding a hard core interaction is not 
easy to handle and more complicated, many body long range interactions have been intro- 
duced (see \13\ I15|) to handle this in the framework of the equilibrium Statistical Mechanics. 

While one could in principle try to follow the approach in [13], we find easier to consider 
a different kinetic model introduced a few years ago in [2], where there is no attractive 
interaction and the Statistical Mechanics instability does not arise. The model consists of 
two species of particles which, for simplicity, have the same mass. To fix the ideas, think of 
them as distinguished just by their color, say red and blue. Their short range interaction is 
modeled by Boltzmann-like collisions which are color blind, while the long range repulsive 
interaction, arising only between particles of different color, is modeled by a Vlasov force with 
a smooth, bounded, finite range potential. We refer to [2], [I], [3], [1] for more information 
on this model. For any t G M"*", let the non negative normalized functions fi{t,x,^), z = 1,2, 
denote the probability densities of finding a particle of species 1 (red) or 2 (blue) in a cell 
of the phase space around the point (x,^) at time t. In this paper we will only consider the 
case x G M, the real line, while the velocity ^ = (u, ^) G with £ R^. The time evolution 
is governed by 

dtfi + vd^h + F{f2)d,h = Qifi, fi + /2), 

(1.1) 

dtf2 + Vd,f2 + F{h)d,f2 = Q(/2, /l + /2). 

Here the Vlasov force F{h) due to the mass distribution h is defined as 

F{h){t,x) = -dJ dyU{\x-y\) f dCh{t,y,0, (1-2) 

where U{r), the interaction potential is non negative, smooth, bounded, with U{r) = for 
r > 1 and 

/ dxU{\x\) = 1. 
Jr 

The collision integral is defined as: 

Q{hi,h2m = I \{^-0-^\{hi{^*)h2{C)-hi{Oh2ia}dC'di^ 

(1.3) 

:= Qgain(^l, ^2) — <5loss(^l, ^2), 

with = {w G I |a;| = 1} and ^'^ related to ^, ^' by the usual elastic collision relations 

= c - • (e - ^% e: = e' + ^[^ • (e - c')]- (i-4) 

For any (3 > 0, the couple {aifii3,a2fii3), with /i^ the spatially homogeneous Maxwellian 
at temperature 

M.(e) = (^)%-^^^/^ (1.5) 
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and aj > 0, is an equilibrium solution and the most general homogeneous equilibrium differs 
from it just for rescaling and centering. However, due to the presence of the Vlasov force, 
non homogeneous Maxwellian equilibria are possible and they are of the form 

h{x,i) = pi{x)^xp{i) (1.6) 

provided that the densities Pi{x) > satisfy the conditions 

logp^ + pU*p^+i = Ci, i = l,2 (1.7) 

Here and in the rest of the paper the label i + 1 means 2 if i = 1 and 1 if i = 2 and the 
convolution product * is defined by a * b{x) = Jj^dya{\x — y\)b{y). The constants Cj have 
the physical meaning of chemical potentials. The phase transition phenomenon we want to 
discuss occurs when Ci = C2, a situation where the equilibrium solutions are symmetric under 
the exchange 1^2. This symmetry is spontaneously broken below the critical temperature. 
Therefore we assume Ci = C2 in the rest of the paper. 

If (3 is suitably small, conditions (jl.7p are satisfied only by constant piS. On the other 
hand if (3 is sufficiently large, non constant pi solving (II. 7p can be constructed. 

To understand the arising of multiple equilibria, let us start by defining the local free 
energy ip{pi, P2) on x M+ (i.e. the free energy density when U is replaced by a (^-function) 
as 

^iPl,P2) = /Ollogpi + p2^0gp2 + (3piP2 (1-8) 

Let p = pi + P2- It is shown in [4J] (in a slightly more general context) that, if /3p < 2 then 
the only stationary points for ip are characterized by pi = p2 = p/2 and they are minimizers 
for ip. On the other hand, if /3p > 2, then there are p^ > p^ > such that 

(,Pi,P2) = {p^,p-) (1.9) 

(and {pi, P2) = {p^iP^) by symmetry) is an absolute minimizer for 99, while (pi,p2) = 
(p/2, p/2) is a local maximizer for ip. When (ai, 02) is chosen equal to {p'^, p~), the equilibrium 
state is interpreted as a pure phase rich of red particles, while the state (01,02) = (p^ , p'^) 
corresponds to a phase rich of blue particles. 

Non homogeneous solutions become relevant when one has to describe a situation of 
phase coexistence. The idea to construct non homogeneous solutions is to observe that at 
low temperature, in order to minimize the local free energy, the system has to be in a pure 
phase at infinity. A situation with phase coexistence should arise when the system is in the 
blue-rich phase at —00 and in the red-rich phase at -|-oo or vice versa. Therefore one looks 
at the minimizers of a suitable free energy functional with above constraints. Indeed, the 
conditions (jl.7p are the Euler-Lagrange equations for the free energy functional 

^{pi,P2) = j^dx{pi\ogpi + p2log/52) + P j^dx j^dyU{\x - y\) Pi{x) p2{y) ■ (1.10) 

on a finite interval / with periodic boundary conditions. Moreover, it can be shown, (see 
[3]), that when (3 is large, with the masses of the two species fixed, there are non constant 
couples [pit P2) giving lower values to the free energy than the constants. Over the real line 
—00 < X < 00, above free energy does not make sense and a more careful definition is required: 
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Let {—i, i) be a bounded interval and set 

^{~e/){pi, P2) = dxip{pi,p2) 
Ji-e,i) 

/3 f 

+ TT / dxdyU{\x - y\)[pi{x) - pi{y)][p2{y) - p2{x)]. 

We define the excess free energy as 

Tipi,p2) :=lim [T^_,^,)ipi,p2)-2iip{p+,p~)]. (1.11) 

Note that, since ip{p'^,p^) = (/?(p^,p+), the functional JT is +00 when p = {pi,p2) does not 
go to pure phases at infinity. Since we are interested at the coexistence of phases, we assume 
that {pi, P2) satisfy one of the conditions lima;^-|-oo Pi{x) = p^ or lim^^^-i-oo P2{x) = p"^ ■ In [3] 
several results are proved about the minimizers for the functional which are summarized 
in the following theorem (see also|5j): 

Theorem 1.1 Let (3p> 2. Then there exists a unique (up to translations) positive minimizer 
(front) for the one- dimensional excess free energy T , defined by in the class of contin- 

uous front functions p = {pi{x), p2{x)) such that lim^^^-i-oo Pi = P^, lim^^±oo /52 = P^ , 
where p^ are defined in lll.9\) . We denote by p = {pi{x), p2{x)) the unique minimizer 
such that pi{x) = p2{—x). pi is monotone increasing and p2 is monotone decreasing and 
p~ < pi{x) < p"*" for any x G M. Moreover, the front p is C°"{M)-smooth and satisfies the 
Euler- Lagrange equations (fi. 7| ); Its derivative p' satisfies the equations 

-p[ + l3-pi{U*p[+^) = Q, i = l,2. (1.12) 

The front p converges to its asymptotic values exponentially fast, in the sense that there is 
a > such that 

\pi{x) - /j'^le"'^' as X ^ ^00, \p2ix) - p'^je"'^' as x ^ ^00. 

Finally, the derivatives of pi of any order vanish at infinity exponentially fast and p' is odd 
in the sense that 

Pii^) = -M-x)- (1-13) 

The aim of this paper is to show that for (3 larger than a certain critical value, the non 
homogeneous equilibria, the front solutions {pi,p2) are dynamically stable with respect to 
the evolution ()1.1|] . while the homogeneous mixed phase is unstable. This will show that a 
complete bifurcation scenario arises in this model for any value of the Knudsen number, i.e. 
the ratio between the mean free path and the range of the interaction potential. 

In the rest of the paper, without loss of generality we fix the asymptotic total density 
p = p^ + p^ = 2. If /? < 1, the only minimizer is the homogeneous equilibrium (mixed 
phase) -Mhom(0 = (P/9(0) P/3(C))- Oii the other hand, if /3 > 1, the pure phases Mred = 
{p'^ Ppii)-, p~ p-liii)) and Mbiue = {p~ 1^/3 iO^ P^ P 1^(0) ^^e constant minimizers and Mhom is a 
maximizer. Moreover, we have the non homogeneous equilibrium 



= {pi{x)pp{i),p2{x)pp{Cl), 



(1.14) 
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with pi{x) + p2{x) ^ 2 as X ^ ±00. 

We assume that the initial datum for the dynamics (jl.ip is a small perturbation of one of 
the above equilibria, denoted generically by M: 

/i(0) = M, + v^5i(0), (1.15) 

with Qi sufficiently small in a sense that will be specified later. 

Let f{t,x,(,) be the solution to (jl.ip and define g{t,x,(,) by setting 

Mt,x,0 = Mi{x,i) + ^Mi{x,i)gi{x,i,t). (1.16) 

The equation for the perturbation g is 

^dt + vd^ + F [Ui+i + VMi+i5i+i) ^1- + ^{x, 0} <?i 

2 

+r(9i,9i) + r(5i,5i+i), 

where we have used the notation: 



Kx,6= / \{i-i')-u:\{M^{i')+M2{i'))di'dL0, 




(1.18) 



We will need the following symmetry condition on the initial data: 

5i(0,x,?;,C) =52(0,-2;,-t;,C). (1.19) 

We note that such a property is preserved by the time evolution. It plays a crucial role in 
removing the obvious orbital instability of our problem related to the translation invariance: 
indeed such an invariance is explicitly broken when condition (jl.l9p is assumed. 

The momentum, kinetic energy and particle number of each species are conserved during 
collisions, while the sum of potential and kinetic energy is conserved along the trajectories. 
Therefore the following quantities are conserved under the evolution (jl.ip : 

• The masses of the perturbation g, 

Mi{g) = [ dx [ d^Mx, - Mi{x, 0), i = h2; 
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The energy of the perturbation g, 

8{g) = [ dx\[ di^ [ih + f2) - (Ml + M2) 



+ / dyU{\x - y\)[pf^ix)pf^{y) - pMAx)pM2iy)] 

JR 

with fi = Mi + \/Migi and 



(1.20) 



Pfix)= / dUix,0, (1-21) 
Jr3 

the spatial density of the distribution /. 
Moreover, by standard arguments it follows that the //-function of the perturbation g, 

H{g) = [ dx [ dd{fi{x,0logfi{x,0+f2{x,0^ogf2{x,0) 
Jr Jr^ l 

- (Ml {x, log Mi(x, + M2{x, log M2(x, 0)] (1-22) 

does not increase in the evolution (jl.ip . 

We notice that, since the spatial domain is M, the total masses, energy and //-function of 
the distribution / are not well defined, while above differences are finite. Since H{g) does not 
increase during the evolution and £{g) and the total masses Mi{g) are constant, any linear 
combination of them, with positive coefficient for the entropy, does not increase. In particular, 
the following non- increasing entropy- energy functional is crucial to study the stability of the 
equilibria: 

n{g) = H{g) + P£{g) - J2M^{g) + 1 + In (^^) ' | . (1.23) 

The factors multiplying the masses are suitably chosen to cancel some linear terms, as will be 
shown in next section. The factor /? in front of the temperature is dictated by the free energy 
minimizing properties of the equilibria. 

In the next sections we shall use the following weighted norm for p G [1, +co]: 

dx [ d^\wiOf{x,CW 
Jr3 

for some positive weight function w. Moreover L^^p denotes the space of the measurable 
functions om M x with || • \\lp,w finite. When notation will be unambiguous, we will also 
omit the index w. Finally 'Vx,v denotes the couple {dx,dy). 

Theorem 1.2 (Stability) Assume (3 > I and M = Mp. Let w{i) = (S + |^p)T, /or S > 
and 7 > 3/2. There are (5 > and S > 0, C > such that, if the initial datura fi{0) = 
Mi + \/Migi(0) satisfies the symmetry condition \1.19\) . the bound 

\\wgm\L^ + Vn{gm < 6 

and II Va;^^,5f(0) 11/^2 < +00, then the initial value problem for (f/./T] ) has a unique global in time 
solution with 



sup \\wgit)\\L^ < C{\\wgiO)\\L^ + ^W(<?(0))}, (1.24) 

0<t<oo 

\\Vx,vgmL^ < e^*||V,,,<7(0)||L2. (1.25) 
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Remark 1.1: Theorem \1.S\ implies the stability of the non homogeneous equilibrium solution 
Mp in the Vlasov-Boltzmann evolution with respect to initial perturbations satisfying 

^.19\) . It turns out that H{g{0)) > for \\wg{0)\\L^ small (see Lemma [2^) . 

Dynamical stability and rate of convergence for the same front solutions has been estab- 
lished [5] for the Vlasov-Fokker-Planck dynamics: the Vlasov force is the same introduced 
here, but the collisions are replaced by a Fokker-Plank operator modeling the contact with a 
reservoir at inverse temperature p. A careful analysis of the macroscopic equation plays an 
important role. 

To prove nonlinear stability, the first major mathematical difficulty we encounter is the 
presence of a large amplitude potential U. To our knowledge, so far there has been no pub- 
lished work on stability result in the presence of a large external field in the Boltzmann theory. 
The main problem is the collapse of Sobolev estimate in higher order energy norms. Indeed, 
even upon taking one x— derivative, the iif^-norm might actually grow in time due to the 
presence of the term dxF{Mi+i)dvgi in ()4.14p . We are thus forced to design a strategy of 
proof based on a weighted L°° formulation without any derivatives and get control of deriva- 
tives only afterwords. Furthermore, unlike the previous linear Fokker-Planck interaction, the 
nonlinear Boltzmann collisions make it difficult to analyze the equation for the projection on 
the hydrodynamical modes in this case. In fact, even a stability of g is difficult to obtain 
directly from the analysis of eq. (jl.l7p . Instead, we avoid a direct study of eq. (jl.l7p and 
make crucial use of the fundamental entropy-energy Tl{g) estimate to obtain a mixed — 
type of stability estimate, based on the spectral gap of the linearized free energy operator 
A (Lemma 12. 2p . We note that the spectral gap estimate relies in an essential way on the 
minimizing properties of the front solution and the symmetry condition is used to control the 
part of the solution in the null space of the operator A. We then bootstrap such a LP' stability 
to a estimate to obtain pointwise stability estimate, by following the curved trajectory 
induced by the force field F(Mj+i -|- ^Mj+i^/j+i). The success of such a strategy is somewhat 
surprising: the perturbation of the field does not need to decay in time, and all the analysis 
over the one dimensional trajectory is carried out in a finite time interval [0, Tq] (Lemma 
14. ip . Yet this is still sufficient for the global in time estimate due to the strong exponential 
time decay from the collision frequency. 

Remark 1.2: The same result also holds when perturbing initially the equilibria Mj-ed o-nd 
-^biue- Indeed the proof is even simpler in this case, because the the analog of the operator 
A has again a spectral gap property, but its null space is trivial. Due to this, the symmetry 
assumption U.19\) is no more needed. 

We next discuss the homogeneous equilibrium Mhom = ilJ'/s, IJ-fs)- 

Theorem 1.3 Assume (3 < I and M = Mhom- There are 6 > and S > 0, C > such that, 
if the initial datum /i(0) = Mj + \/Migi{0) satisfies the bound 

Wwgoh^ + \/W(5(0)) < 5 

and \\Vx,vg{0)\\L'^ < +oo, then the initial value problem for Ji. j7| j has a unique global in time 
solution satisfying the estimates ll.24\) and U.25\) . 

Remark 1.3: This implies the stability of the mixed phase for [3 < 1. The proof is again a 
simpler version of the one for Theorem \L2\ which will be omitted. In this case the analog of 
operator A has a spectral gap property (because (3 < 1) and a trivial null space, so that we do 
not need to require the symmetry property U.19\) . 
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We have already noticed that, when (3 > 1, the couple (1, 1) is a local maximizer for the 
free energy. In next theorem we establish the dynamical instability of Mhom- 

Theorem 1.4 (Instability) Assume /? > 1. There exist constants ko > 0, 9 > , C > 0, 

c > and a family of initial "^—periodic data //(O) = ^i3+yJJjLpgf{^) > 0, with g^{0) satisfying 
Kl^) and 

II V.,./ (0) 11^2 + (0)11 LOO <C5, 
for 5 sufficiently small, hut the solution g^{t) to li. j|) satisfies 

sup \\wg\t)\\L^>c sup ||/(t) 11^,2 > 06* > 0. 

0<t<T^ 0<t<T^ 

Here the escape time is 

A is the eigenvalue with the largest real part for the linearized Vlasov-Boltzmann system con- 
structed in Theorem \5.S\ with Re A > 0. 

Remark 1.4: Note that ^ oo as 6 ^ 0. We also observe that the critical value (3 = 1 
escapes our analysis, since it is based on some strict inequalities which collapse for the critical 
(3. Furthermore, the growing mode which we construct satisfies the symmetry condition il.l9\) . 
Hence the instability is not a consequence of the absence of such a symmetry. 

To prove instability of the homogeneous state, we encounter a second major difficulty. 
Even though such a homogeneous equilibrium is not a local minimizer for the free-energy 
(11. lip , it is not clear at all if such a property leads to a dynamical instability along the time 
evolution. As a matter of fact, the presence of (possibly strong) stabilizing collision in the 
velocity space might damp out the instability. It is thus fundamental to understand such a 
stabilizing collision effect. Unfortunately, a direct linear stability analysis in the presence of 
the collision effect is too complicated to draw any conclusion. We therefore developed a new 
perturbation argument to establish the linear instability around such a homogeneous steady 
state. We first observe that in the absence of the collision effect, the homogeneous state is 
indeed dynamically unstable, by an explicit analysis similar to the Penrose criterion in plasma 
physics [14J. It then follows, via a perturbation argument, that for 'weak' collision effects, 
such an instability should persist. The key is to show that this is even true for arbitrarily 
'strong' collision effect. We use an argument of contradiction and a method of continuation. 
In fact, if instability should fail at some level of the collision effect, then a neutral mode must 
occur. The interaction between the Vlasov force and the collision effect forces any neutral 
mode to behave like a multiple of Mhom with a particular dispersion relation never satisfied 
for our model (Theorem 15. 2p . To bootstrap such a linear instability into a nonlinear one is 
delicate due to severe nonlinearity, and we follow the program developed by Strauss and the 
second author over the years [lOl E] • 

We remark that the instability of the homogeneous equilibrium for VFP model is still open 
because the techniques used in the present paper cannot apply due to the unboundedness of 
the Fokker-Plank operator. On the other hand, in this paper we cannot prove the convergence 
to the stable equilibrium, while it was possible in the VFP model even to compute the rate. 

The paper is organized as follows: In Section 2 we use the energy-entropy to derive a mixed 

— estimate. In Section 3, we establish some lemmas on the characteristics curves for the 
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equations (jl.ip . In Section 4 we establish the nonhnear stabihty in L°° norm. In Section 5 we 
construct a growing mode for the hnearized problem around the homogeneous equilibrium and 
finally, in Section 6 we show that such a linear growing mode leads to nonlinear instability. 



2 Entropy-Energy Estimate. 

In this section we use the conservation of energy and masses, and the entropy inequality to 
obtain a priori estimates on the deviation of the solution from equilibrium. A crucial role is 
played by the quadratic approximation of the excess free energy functional T. We need a few 
definitions and some notation. 

Let u = {ui,U2) be a couple of functions in L2(M) and denote {u, v) = X^i=i('Ui, t'i)- Given 
a couple of densities p = {pi{x), p2{x)), we define the operator A by setting 



{u,Au) ^ y dxui{x){Au)i{x) = 2^-^( P + ^'")L=o 



In particular, when p = {pi{x), p2{x)), the action of the operator ^4 on u is 

(Au)i = — + pU*U2, {Au)2 = — +PU*ui. (2.1) 

Pi P2 

Hence 

{u,Au) = ^ V / dx^^^+^ [ dx f dyU{\x -y\)ui{x)u2{y). 
2^ 7m Pi{x) Jr Jr 

Due to the minimizing properties of p, this quadratic form is non negative. Moreover, by 
(|1.12|) and (j2.ip . we see that Ap' = 0, which shows that p' is in the null space of A. Indeed, 
one can show (see [5j and references quoted therein) the following 

Lemma 2.1 Suppose (3 > 1 and p = {pi,p2). Then there exist 5o > such that 

{u,Au)>6o{{I-V)u,{I-V)u), 

where V is the projector on Null^.' 

Null A = {u£ L2(IR) X l2(R) |n = cp ,c G M}. 

// either p = (p^, p^), or p = (1, 1) with (3 <\, we have 

(n, Au) > 6o{u, u), 

and the null space of A reduces to {0}. 



We prove the following lemma which plays a crucial role in the proof of the stability of 
Mp. Remind that we have adopted the notation ^ = {v,(), with ^ G and G R for the 
velocity. 
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Lemma 2.2 Let M = Mp. Let fi{t,x,v,C) = f2{t, —x, —vX) ctnd that \\g\\L°° ^ ^ for some 
small 6. Then there exists C > and k > 0, such that 



fi{t) - M,|l|^^(i)_M,|>.M,:| < n{gm. 

Remark 2.1: Note that, when dealing with Mred; -^biue cind Mhom; the functions pi and p2 
have to be replaced by the constant values (p'^,p~), {p~ , p'^) and (1, 1) respectively. With this 
modifications the lemma still holds. 

Proof. Remember the notation pf{t,x) = d^f{t,x,(^). We may construct solutions (see 
[8]) such that 

n{g) < n{gm. 

We expand TCig) and use ()1.7p to cancel the linear part of the expansion, which takes the 
form 

2 

- V / / dC{Ci + l + ln(^)i}{/, - Mi} + 

V / / - Mi) + {In Mi + l}{fi - Mi) 

~lJR JR^ [ ^ 

+/?{?7*p,+i}(/i-Mi)l. 



Indeed, since In Mj = In(^) 2 — by (jl.7p . the above quantity is zero by construction. 
Therefore, we turn to the second order expansion of TCig). For some fi between Mi and fi, 



no) = EI di ^^''^'\^^'> dxdi (2.2) 

^ IR Jr3 2fi 



i=l 



+(3 dx dy{pf^{t,x) - pi{x))U{\x - y\){pf^{t,x) - p2ix)). 

JR JR 

For some small number k to be determined, we introduce the indicator functions xf 
hMt)-Ah\<KAh and xf = l|/,{t)-Af,|>KM, and split the first term into 



(/,(t)-M,)2 ^ (■ (■ (fi{t)-Mi^^ 



dx I d c^'{-' xf + I dx I dg -'-; - x< 

Jr3 2fi Jr Jr3 2ft 

We first estimate ^'^^^ ^ — = — ^— in the case of \fi{t) — Mj| > «;Mj. Notice that either fi > 

2fi 

(1 + k)M„ or /i < (1 - K)Mi. If fi>{l + K)Mi, fi{t) < fi, and we have 
\fi{t) - M,\ ^ \fi{t) - M,\ _ ^ M,^^ 1 _ K 



fi{t) - fi fi ~ l + K l + K 
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In the second case fi < (1 — K.)Mi, fi{t) < Mi and 

|/.(t)-M.|^LM!)-M^ = l_A>i_{i_.} = .> ^ 



fi{t) ~ Mi Mi- ' ' l + K 

Combining these two cases and noticing fi < (1 + K,)Mi for \fi{t) — Mi\ < KMi, we conclude 



'^fi Jr Jr^ 2fi 

)-^0\ ,< , ^ 

2(1 + K)Mi 2(1 + k) 

(2.3) 

1 11,1,^0^ 1^ 



2(1 + k) 



f [ dx [ d^g^ixf + ,777^ I dx f d^\m-M,\xf 
J Jm. Jr3 2(1 + Kj Jus 



^ 2( 



/ / dxn^+ 1^ j dx j di\m-Mi\x 

1 + ^) J JR 2(1 + k) 7k 7ir3 



> 

i 1 



where we have set ni{t,x)yjjr^ = P[(/i(t) — Mj)xf], and P denotes the L| -projection on 



We now spht the potential contribution in (12. 3p to get 

(5 [ dx [ dy [ di f dv^/WigixfUilx - y\)^/M'2g2X2 



+(3 [ dx [ dy [ di f dr] ^/WmxfUilx - y\)^/M~2g2X2 

JR Jr JR'i JR3 

+P I dx I dy I dU drj ^/M[giXiU{\x - y\)^292X2 
JR JR JR^ JR3 

+P [ dx [ dy [ di I dTj ^/M'mxiUi\x-y\)^/M'2g2X2- 

JR JR JR^ JR3 

From our assumption H^Hloo < 6, the last three terms are controlled by 

cJWgiU^ + II52IIL-) dx d^\fi{t) - Mi\xf 

< CpSy^ [ dx [ dC\fi{t) - Mi\l\f^(^t)-M,\>.M,- 

■^^ JR 

which is bounded by the second term in (j2.3p for 5 « k. Since / d^y/Mlgi = ^/pini, the first 
term can be written as 



P dx dy ni{t,x)y/pi{x)U{\x - y\)n2{t,y)^/ p2{y) 

JR JR 
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We now combine it with the first term in (j2.3p to get 

1 



2(1 + k) 



dx{ni + n^) 



+(3 dx dynl{x)^/pl{x)U{\x - y\)n2{y)\/ P2{y)dxdy 



= {n^,An^)-—— — - / dx{nl+nl). 
2(1 + k) Jji 

Since ni(x) = n2{—x) by our symmetry assumption and from p'i{x) = —p2{—x), it follows 
that (ni, is orthogonal to the null space of A. From the spectral inequality for the operator 
A, since Pi{x) > p~ , the above is bounded from below by 

^~^°-2(lT^) + 

2 

> ||P(/i(i) - Mi)lif^^t)_M,\<KMi\\l2, 

1=1 

provided that we choose k sufficiently small. Then the lemma follows by collecting the terms 
with S « K. ■ 

3 1-Dimensional Characteristics. 

We define the characteristics curves [Xi{s;t,x,v),Vi{s;t,x,v)] for ( I1.17P passing through 
{t, x, v) at s = t, such that 

dXi{s;t,x,v) irr ^ \ /on 
— = Vi{s;t,x,v), (3.1) 

dV(s' t x v) r 

y = -d,U* di{Mi+i + yM;y^g,+^) = -d^^{Xi{s-t,x,v)). 

ds Jk3 

We also define the unperturbed characteristics [Xf{s;t,x,v),V^{s]t,x,v)] passing through 
{t, x, v) at s = t, such that 

dX'^{s;t,x,v) _ 



ds 



Vl'{s;t,x,v), (3.2) 



dV^{s;t,x, 



V 



ds 

dXi{s] t, X, V 



Our main goal is to study the zero set of 



dv 



dX'^is; t, X, v) 

Lemma 3.1 For any {t, x, v) with v ^ 0, the set of{s £ M : — — - — = 0} is countable. 



dv 



Proof. From the particle energy conservation for (j3.2p : 



hv^is; t, X, v)\^ - MX^s; t, x, v)) = \\v\^ - ^x). (3.3) 
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Taking derivative with respect to v yields 

Vi [S,t,X,V) (po - V. 

dX'!'(so;t,x,v) 

Assume v ^ 0. If there is sq such that — ^ — ^ — = 0, then necessarily 

dv 

Vt [so;t,x,v) — - — — / 0, 

because, for such an sq it reduces to v by previous equation and we assumed v 0. Hence 

dVP{so;t,x,v) d dX^{s;t,x,v) 



dv ds dv 



=^0 



dX'^is; t, X, v) 

and therefore — — ^ — 7^ in a neighborhood of s = sq. This implies that for v ^ 0, 

dv 

, (9X?(s; t,x,v) . 
the set -^s : — = 0[ is countable. ■ 

dv 

Lemma 3.2 Fix Tq > and N > 0. Let \v\ < N. 

1. For any e > 0, there exists sufficiently large so that, if\x\ > Lgr, then for < s < t—e, 

dX^{s;t,x,v) £ Q 
dv 2 

2. For any r] > 0, there exist P finite points \xi.\ < (1 < k < P) and corresponding 
open sets 

= [J {a-n < s < b„} X {co <v < do} 

(n,o)eIk 

with the property 

\[0,To]x{\v\<N}nOl\<V, 
so that there exists m > and, for any \x\ < L^, there exists I G {!,..., P} 



dXf{s;t,x,v) 



dv 



> m > 0. 



for {s,v) G Oa;r 

Proof. For any e > 0, from l^^i) and ([33] ) 



\VP{s;t,x,v)\ < +2^/1100^- <N + C, 
\Xf{s;t,x,v) -x\ < To{N + C}. 

By choosing (depending on Tq and N) large enough, for \x\ > L^, \v\ < N, and < s < Tq, 

\Xfis;t,x,v)\>^. (3.4) 
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Prom (|3.2p . we have 

dX^{s;t,x,v) 
ds"^ dv 

and we deduce that for |s| < Tq, 



-dxx<Po{Xf{s;t,x,v)) 



dX^{s;t,x,' 
dv 



dX^is;t,x,v) 



dv 

By the Taylor expansion for s, we get 

dX^{s;t,x,v) dX^is;t,x,v) 



<Cto. 



dv 



dv 



^ d dXfis;t,x,v) 

s=t ds dv 



s=t 



£ dXf{s;t,x,v) 
2 ds^ dv 

{s-tf d2 dX^{s;t,x,v) 



{s-t) + 



ds^ 



dv 



(3.5) 



(3.6) 



for some t — Tq < s < t. Since the densities pi tend to the their asymptotic values at infinity, 
limy^oo dxxfpoiu) = 0. Therefore, using again p.5p with s = s and ()3.6p . we have 



dXfis;t,x,v) 
dv 



< {s - t)(l - (t - s)Cto sup \dxx4>o{y)\) 

\y\>¥ 

s - t 



< 



<0, 



by choosing sufficiently large. Part (1) thus follows. 

To prove part (2), for < L^, introduce the zero set of g^*'^'""* as 

^ r ^ ^ , , dXf(s;t,x,v) 
Zx = {t-To<s< To, \v\<N: ' = 0}. 



dv 



Then from the Pubini Theorem and Lemma |3.H 



N 



-N [Jt-To 



axO(s;t,a:,v) ds } dv = 0. 



=0} 



Therefore, there exists an open set Qx such that Zx C with \Qx\ < ^- Clearly g^*''^'""* 7^ 
over the compact set [i — Tq, t] x {|f | < N} n W^. By the continuity in s and f , there exists 
rux > such that over [t - Tq, t] x {|i;| < N} D Qx- 



dXf{s;t,x,v) 



dv 



> Arux > 0. 



Purthermore, from the continuity in x, we have an open set of (x — Ax,x + A^,.) such that 

dXf{s;t, x',v) 



dv 



> 2mx > 



forallx'e {x-Ax,x+Ax), {s,v) e [t-To,t\x {\v\ < N}nQ.x- Such (j;-Aa;, j;+Aa;) forms an 
open covering for |x| < L^, hence there is a finite subcovering {(x^— A^, x^+A^), k = 1, . . . , P} 
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for |x| < Lg. For any |x| < L^, there exists ^Ixi such that x G {xi — Ai,xi + A^) and for 

[t-To,t] X {\v\ <N}nn 



c 

Xl 1 



dXf{s;t,x,v) 



dv 



> 2m = min 2mk > 0. 

1<A:<P 



We finally choose an (finite) open covering of [t — Tq, t] x < N} n of the form 
= Un oi^n < s < bn} X {co < V < do} with |a„ — bn\ + |co — do\ sufficiently small so that 
over Ox^^ 

dXf{s;t,x,v) 



dv 



> m> 



Lemma 3.3 Fix To > and N > 0. Let \v\ < N. For any e > 0, recall L^, Ox^, 1 < k < P 
constructed in Lemma \3.2\ . There exists 6 > such that if WqHl'^ < 

1- If \x\ > then for < s < t — e, 

dXi{s;t,x,v) e 
dv ^ ~2' 

2. For any \x\ < L^, there exists I such that x £ {xi — Ai, xi + A^) and for all (s, v) £ Oxi 



dXi{s; t, X, v) 



dv 



m 

>->0. 



Proof. Denote the solution operator of (|3.2p and (j3.ip by Go and G. By the Duhamel 
principle, we have 

Xi{s;t,x,v)\ ^ ^Go(s-t)fx 
Vi{s;t,x,v) J \v 





Go{s-t) I " (ny) 

dxU*{f ./M-^gi+id^'}{T) I ■ ^ ■ ^ 



Notice that 



d_ 

dv 



dxU*{ y^M;^ig,+ide}{T) 



<C\\g 



dxxU * { J ^/M;^g,+ide} 

dXi{T;t,x,v) 



, dXi{T;t,x,v) 



dv 



dv 



It thus follows from taking derivative of v in (j3.7p and by the Gronwall's lemma that < s < 
To 



dXi{s; t, X, v) 



dv 



< e 



CTo 



We now use p.7p again to get 

dXi{s;t,x,v) dXf{s;t,x,v) 



< C'toI|9||l2- 



dv dv 
Hence, we deduce our lemma by choosing \\g\\L2 sufficiently small. 
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4 Weighted L°° Stability. 

In this section we use the entropy-energy bound and the estimates on the characteristics to 
show that the perturbation g of the non homogeneous equihbrium Mp, is arbitrarily small 
at any positive time in a suitable weighted L^o norm, provided that it is initially sufficiently 
small, thus showing the stability of the non homogeneous equilibrium. We use the weight 
function w(^) = (S + with S a positive constant to be chosen later and 7 > |- 

Lemma 4.1 Let h = wg. There exist Tq > and 5 > such that, if < then 

Proof. We first write the equation for h = wg from (jl.l7|) : 

(dt + vd, + F{M,+i + 01^^^ + u{x, o) hi = 

+F ( Tm^^I vh. + wT(^,^\+wv(^,^\, (4.1) 



W J \W W J \w w ^ 

where K'M-) = wK'^ (-]. We note that for S > 1 

He) ~ WTWW - ' — [s + iep]7 — <c.[i + \^- ei r. (4.2) 

For any (t,x,^), integrating along its backward trajectory (j3.ip 

[Xi{s),Vi{s)] = [Xi{s]t,x,v),Vi{s]t,x,v)] 
we can express hi{t,x,S,) as 

hi{0, Xi{0; t, X, v), Vi{0- 1, X, v),C) + 

ei:Mr)dr^j,^M,+i + y/M;^,g,+,)-h,}is,X,{s),V,{s),C)ds 
w 

+ [ e/>»M^^{F(yii^5i+i)y,(s)v^}(s,Xi(s),yi(s),C)ds 



+ E fef>^^^^''^ I Kii^hA {s,X,{s),Vi{s)X)ds 

+ f e/>»M^^{F(yiV^5i+i)^;/i,}(s,X,(s),yi(s),C) 
Jo 

+ re//-M^-.;[r(^,^) +r(^,^)l(.,X.(.),y.(.),CM.. (4.3) 
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We have set ^'^(t) = z^(Vi(r), C) ^ '^o > . Fix a small constant e > 0. We can choose S 
large so that |^| < e. Since ||-F(Mj+i + ^Mj_|_i(7i+i)||ioo < C if is small, the second 

term in ()4.3p is bounded by 

Cee 2" sup {e~||/i(s)||Loo}. 

0<s<To 

For the third term in (14. 3p . we split 

= 5i+ll|/,(t)-M,|>KMi + 5'j+ll|/,(t)-Mi|>KMr 

Since C/ is smooth, 



\F{y/M;^igi+i)\\i 



< C{|| Y^Mj+i5rj+il|j^(i)_^.|<^;(^J|i2 + II V Mj4.i5rj+il|j.(f)_j\^.|>^j\^. II^,! }, 



by Lemma [27 



(4.4) 



For the fifth term, we note that, since for hard spheres fi(s) > 1^0(1 + ICI + l^(s)('5)|! it 

*d5^lMe//(-.(-)/2)rf- < 1. Therefore, 



\Vi(s)\ 1 
follows that — -— < — . Moreover 



<Ce"''«* sup {eT||/j(s)||^oc}2. 

0<s<To 



For the last term in (j4.3p . by Lemma 10 of [7J, it follows 

hi /ij+i 



wT 



hi hi 
w ' w 



(6 



+ 



wT 



w w 

We therefore get the bound for the last term by 
ft 



(0 



<cumh\\i. 



,I>^^^''^Uiis)\\h{s)\\l^ds 







<C{ sup e^||/i(s)||Loc}2 / e^>»(")'^^z^i(s)e-'^°^ds. 
o<s<ro Jo 



Note that ^[e//^«W'^^] = e^^'i^)^'' Ui{s) . Integrating by parts yields 

< C(l + t)e-'^«*. 



We shall mainly concentrate on the fourth term in ()4.3p . Let k*'-'(^,^') be the corre- 
sponding kernel associated with in (j4.ip . We now use (j4.3p for hj(s, Xi{s),S,') again to 
evaluate 
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Denote 

[Xj{si),Vj{si)] = [Xj{si;Xi{s;t,x,v),v'), Vj{si; Xi{s;t, x,v),v')]. 
We can bound the fourth term in ()4.3p by the sum on j of 



j:u,(r)dr+fy,(r)dr f (^V,{s) , ( , Olhj {0, X, (0) ,V, (0) , (' M' ds 



t 



Jo Jsi J R3 

w' 

x{F{Mj+i + ^Mj+igj+i) — hj}{si,Xj{si),Vj{si),C')dedsdsi 

Jo Jsi J R3 

x{F(VM~^5j+i)i;V^}(si,X,(si),F,(si),C')<c?5 
X ki'=(y, (si ) , C' , e") I ^fc (^1 ,Xj{si), e") I di'd^dsdsi 

Jo Jsi J 



x{F{^Mj+igj+i)vhj}{suXj{si),Vj{si),C')d^'dsdsi 

+ [ r el>dr)dr+j:^ u,ir)dr f | kM (y. , ^, | (4.5) 



'0 Jsi 



xw 

'WW/ \w w 



We will make an extended use of Lemma 7 of which we report here for reader's 
convenience: For hard spheres, the usual Grad estimates imply: 



lit t/|2 1 iiei^-ie'i^P 



\^i,M,a\<m-c\ + \^-cr'}e ^u^. (4.6) 

Lemma 4.2 (Lemma 7 of [Tj) There are e > and C > such that 

/ ^^'^^1^ - ^'1 + 1^ - ^r^}e-^'^-^'' -^^^ < (4.7) 
By Lemma 14.21 we obtain the crucial estimate 

K^\{^,em' <-r^ (4.8) 
1 + 141 



uniformly in E. Since > vq, by taking norm for h and (j4.8p . we bound the first term 
in (14. 5p by Cte~'^''*||/io||L°° , and the second term by 

eCe-^' sup {e^'||/i(s)||Loo}. 

0<s<To 
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By (|4.8p . the third term is bounded by CTC{g{0)) as in (j4.4p . and the last two nonhnear terms 
are bounded by 

C{l + t}e-^«*{ sup e^^||/i(s)||^oop. 
o<s<ro 

We now concentrate on the fourth term in (j4.5p . which will be estimated along the same 
lines of the proof of Theorem 20 in [7] . 

CASE 1: For |^| > Nto, we know that from jM]) 



\Vi{s;t,x,v) -v\< \s-t\C < CTo, 
By Lemma 14.21 and (|4.2p . for Ntq large 



< ^ < ^ 



l + \Vi{s)\ + \C\ - N-CTo' 
we therefore can find an upper bound for the fourth term in (j4.5p by {N >> Tq) 



N 



JO 



^ o<s<ro 



Ce 2 ^ 
< — — — sup e 2 *||/i(s)||2,oo . 



CASE 2: For |^| < N, |^'| > 2iV, or |^'| < 2N, |^"| > 3N. Notice that we have either 
W -^\>Nov If -C'\> N. This implies that 

\v' — Vi{s;t, x,v)\ > \v' — v\ — \v — Vi{s;t, x,v)\ 

> \v' -v\- CTo, 

\v" - Vj{si;Xi{s;t,x,v),v')\ > \v" - v'\ - \v' - Vj{si; Xi{s;t, x,v),v')\ 

> \v"-v'\-CTo. 

Therefore, either one of the following are valid correspondingly for some o" > 0: 

\ki;,HVj{si),C',n\ < CToe~f^'|ki'=(y,(.i),C',nx 

e|{|V,(si)-^"|2+|C"-CT}|. 

From Lemma l4.2[ 



+ I |ki;^(y,-(.i),C',e")et^l^^(^^)-^"l'+l^"-^'l'>K<+oo. (4.9) 
We use this bound to combine the cases of |f — f > or |f — f | > as: 

ft fSl 




+ 

/O Jo \J\(\<N,\C\>2N, JW\<2N,\("\>3N 
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We first integrate ^' for the first integral and apply (j4.8p to integrate k^^*^ over ^" . We then 
integrate ^" for the second integral and apply ( 14. 8p to integrate k^-* over We thus find an 
upper bound 



c[ Tsup/ \k]i{V,is),C,a\de 

Jo Jo $. J\(\<N,\^'\>2N, 

+ sup / \k^{V,{s,),C',em" (4.10) 

£' J|f'|<2Ar,|f"|>3Af 



Jo Jo 

< C^e-i^'e-^* sup {e^'\\h{s)\\L^}. 

0<s<t 

CASE 3. 1^1 < N, \^'\ < 2iV, |^"| < 3A^. This is the last remaining case because if 
IC'I > 2N, it is included in Case 2; while if |^"| > 3A^, either |^'| < 2N or |^'| > 2N are also 
included in Case 2. We now can bound the second term in (j4.5p by 

c f [ /%^'^''(*-^i)|kJ^^(y,(s),c,eOki'(vs(5i),C',C'')/ifc(si,^.(^i),C'')l 

JO Jb Jo 

where B = {|^'| < 2N, \^"\ < 3N}. We notice that k^^'(^,^') has a possible integrable 
singularity of the type We can choose k^(^,^') smooth with compact support such 

that 

\p\<3NJ\a<3N 



sup / |k'^^(p,^')_kji.(p,e')|d^'<^. (4.11) 

A<3N J\^'\<3N 

Split k^(y,(s),C,e')ki'(^,(si),C',e") into 



We then integrate the first term above in ^" and the second term above in By (j4.8p . we 
can use such an approximation (j4.1ip to bound the si,s integration by 

Ce 2 '■ 

— — sup {e^'\\his)\\L^} 

0<s<t 

x| sup /|kj;^(y,(si),c',e")K+ sup [\yj'ms)x,a\de}\ (4.12) 

[\^'\<2N J \i\<2N J J 

+c f [ re-'^"(*-^^)k^(y,(s),c,e')4(^j(^i),e")|/i.-(s,^.-(si),c',e")|. 



JB Js 



The first term above is further bounded by — supo<5<^ 

Fix e > 0. We use now Lemma 13.31 for the last main contribution in (j4.12p for 
which we separate two cases \Xj(s;t, x,v)\ > and \Xj{s;t,x,v)\ < where is given in 
Lemma 
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In the case \Xj{s;t, x,v)\ > L^, we bound it by 

Cn f I re-^»(*-^^)|/ifc(s,X,(5i),e")|l|x,(s,)|>L.d^de"<dsi 

Jo JB Jo 



10 JB JO 

< Cn{ + / }. 

Jo Jt-e 



The second integral is bounded by Cee 2 *supo<s<t{e 2 *||/i(s)||loo}. In the first integral, 
since s <t — £,hy Lemma 13.31 we can make a change of variable 

y = Xj{si) = X.j{si-Xi{s;t,x,v),v') (4.13) 
because |^| > |. We observe since < C, that from (I3.1jl 

\V -V,{t)\ < J^\\d^nL^dT<To\\d.,<P\\L^, 

\y-X^is)\ < r\Vj{T)\dT<Toi\v'\+To\\d^(P\\L^)<CTo,N 

J SI 

for < 2N. By first integrating over and using the change of variable (j4.13p 

^'''^'-''^\hk{sux,isi),a\'i-\x,(s,)\>LA'^dededs, 

r V / r e-'''^'-'^^\hk{si,y,e)\dsd^"dyds, 

Jo J\y-X,{s)\<CTo,N J\i;"\<3N Jo 

< — sup [ [ \hk{si,y,i")\dsidi"dy 

e 0<si<To J\y~-X,{s)\<CT^^M JW'\<^N 



c 

< — 

e 



t—e r- f-s 
/ / 

Jb Jo 

t—e r f fs 



+ 

I /fc (t) - A4 1 > kM, J I /fe (t) - Mfe I < kMj 

< CT,,NAn9io)) + Vn(m)}- 

We have used the fact = — ^"^^7^^, (which is bounded by fk — for |,^"| < 'iN), and 
applied Lemma |2. 21 

For |Xj(s; t, X, f)| < L^, for any r/ > 0, we again employ Lemma [3^ to find such that 



+ 



[ " / [ \-'''^'-'''>\hk{suX^{si),e)\l\x^(^,)\^Ljsded^'dsi 
Jo Jb Jo 

/ / / lose-^''(*-^^)|/ifc(si,X,(5i),e")|l|x,(s)|<L,c^sde"de'dsi 

JO J B Jo 

[ ° / / \o.e-'^«(*-^^)|/ife(5i,X,(5i),e'0|l|x,(s)|<L,d^^ieX'd5i. 
JO Jb Jo 



Since | [0, Tq] x [-N, N] nO^ \ < t], the first part is bounded by 



CTo,N,er]e 2 * sup {e 2 *||/i(s)||loo}. 

0<s<t 
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The second part is bounded by 



CTo,N,e r r I lo \hk{si,Xj{si),e')dsds,dede'. 

Jo Jo Jb 

Since 1 9Xj{s-i,xj^s^,t,x,v),v ) | ^ ■fji^^2 on from Lemma 13.31 we can make a (local) change of 
variable y = Xj(si) = Xj{si; Xi{s;t,x,v),v') to get 

CTo,N,e / / / ^O^Jhk{si,Xj{si),(,")dsdsid(,'d(," 
Jo Jo Jb ' 

Jo 



X 

J\y-X,(si)\<CTo,N 



hk{si,y,S!')dsdsidydi" 

'|g"|<3Af 

- I 

J\fk{t)-Mk\<KMk J\fk{t)-Mk\>KMk 
< CTo,N,e,r,mgm + ^^(^(O))}. 

Collecting terms, we conclude 

sup e'?*||Ms)l|L- <C(l+To)||/i(0)|U=o 

0<s<To 



+{^ + CN,Toe + Cn,To,sV} sup {e^'\\h{s)\\L^} 

o<s<ro 



+C{ sup e^mhis)\\L^}^ + CTo,N,eWn9m. 

0<s<To 

Assume supo<s<7|, is sufficiently small. We first choose Tq sufficiently large so that 

2C(l + ro)e-^^« < i 

then N sufficiently large, then e sufficiently small, finally rj small to conclude our lemma. ■ 
Proof, of Theorem 11.21 Assume supo<«oc ||/i(t)||Loo is small. We first establish (jl.24p . 
Choose any n = 0, 1, 2, 3, . . . and apply Lemma [4T] repeatedly to get 

1, 



\\hinTo)\\L^ < -\\hi{n - 1}To)\\l^ + CToVngm 

< ^\\h{{n-2}To)\\L^ + ^CToVn{m) + CToVn^ 

< ... 

< ^\\ho\\L^ + Cro VW(9(0)){1 + \ + \ + ---} 



< -\\ho\\L^ + 2CToVnm)- 

For any t, we can find n such that uTq < t < {n + l}Tb, and from L°° estimate from [0,To], 
we conclude (|1.24|) by 



ii/i(t)iiL- < CTo\\HnTo)\\ < c{\\ho\\L^ + vnm)}- 
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To prove (|1.25p . we take x and v derivatives to get 

{dt + vd^ + F{Mi+i + y^m^gi+i)d, + u{0}d,gi - K^'^d^gj 
= -d^F{Mi+i + y^Mi+igi+i)d^gi + /3d^F{y^M~:^9i+i)vd^VMi + 

+d^{F{y^g,)vgi} + dAr{gi,gi) + T{gi,gj)}; (4.14) 

{dt + vd^ + F{Mj + ^/M'jg,)d, + u{0}dvgi - d^K'^^gj} + {dM0}9i 
= -d^9i + Pd^F{y^jgj)vy^, + f3d^F{y^jgj)vd^y^i 

+F{y^g,)dAvg,} + dA^{9i,9i) + ^{9i,9j)} (4.15) 

where K^'^ has similar property as Ki in [6j (see Lemma 2.2 in [6J, p. 1109. In particular, 
\\dv{K^'^ gj}dvgi\\ii < ^\\dvg\\1 + C\\g\\'^j^2 so that a positive dissipation for d^gi occur for small 
||/i||2,oo in (|4.15p . Notice that L = — K > 0. We take inner product with dxgi and d^gi 
respectively, following the procedures in ^ to get: 

f^lWd.gWh < C{||a,F(M,+i)||ioo + \\h\\L^}\\V,,,g\\l, + C\\g\\l,. 
JtlW^^aWh + \\\d.9\\l < C\\d,9\\h + C\\9\\h- (4.16) 

Hence (I1.25P follows from the Gronwall Lemma since supo<t<oo ll^(*)ll-L°° is bounded by (jl.24p . 
With such an estimate, we obtain the uniqueness by taking estimate for the difference for 
(|1.17p because the most difficult term F(Mj+i + ^/Mi^gi^i)dvgi can be handled. ■ 



5 Linear Instability: Growing Mode. 

In this section we study the linearization of the equation (jl.ip around the homogeneous 
equilibrium Mhom = (/^/3)/"/3)- In the sequel we omit the index (3 for sake of shortness: 

/u = /i/3 = (^)^/^e"'^~. When M is replaced by Mhom = (/^, A*) in (11.17p . we get the 
following linearized Vlasov-Boltzmann system: 

dtg + j:g = 0, (5.1) 

where 5 = (51,92), 

{Cg)i = vd^gi - f3F{^gi+i)vy/Jl - Ug 

and ^ 

Lig = — {Q{^/Jigu 2fi) + Q(^, ^(c/i + 5-2)) 

We seek an exponential growing mode for such a system when /? > 1. To this end, we consider 
a family of systems 

dtg + £"5 = 0, (5.2) 

(£" g)i = vd^gi - (5F{^gi+i)v^ - aUg 

and show that there is a growing mode for all a > O.We shall first seek a growing mode of 
the form 51 = 52 = 9, so that the system (|5.2p reduces to the single equation 



{dt + vdx}g - (3F{y/Jlg)v^/Jl = aLg. 



(5.3) 
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with Lg = (Qi^/Ji-g, lA + \^9^ ■ We seek a growing mode periodic in x, so we assume 

periodic dependence on space and exponential in time: g{t,x,S,) = e'^^ e^^^ q{£) . From the 
definition of F (see (|1.2p ). (j5.3p becomes: 

{A + ivk}q - (3k\U{k) | j qy/Jld^^ Vy/Ji = aLq. (5.4) 

Equivalently, q((^) is an eigenfunction for the operator T" with eigenvalue A : 

(T"<7)(0 = ivkqiO - (ik\U{k) I j q{i)^di^ - aLq{i). (5.5) 

Lemma 5.1 Let (3 > 1. There exists sufficiently small a > such that there is an eigen- 
function q{^) to T" with ReA > 0. 

Proof. We first study the eigenvalue problem for the unperturbed operator for a = : 

{X + ivk)q- pkiU{k)!^j q^d(^v^ = 0. (5.6) 

This is similar to the Penrose dispersion relation for the unperturbed Vlasov-Poisson system 
(|14j) for a collisionless plasma. From (j5.6p . we obtain: 



q 



PkU{k)iy^^q^d^^v^ 



(A + ivk) 



Normalizing / q^/Jld^ = 1, we deduce that 

JR3 



Jr3 (A + ivk) 

Multiply and divide by (A — ivk), take the imaginary part and then integrate on ^. By 
consistency we must have: 

Define F(\, k) = P [ ^^^(^)feVO ^g^ Clearly F(0, k) = (3Uik). Since C/(0) = 1 and /? > 1, 

7r3 A^ + kH^ 
there is feo sufficiently small so that 



F(0,A;o) =/5C/(A;o) > 1. (5.7) 

Moreover, lim^-^oo -^(A, feo) = for any ko ^ 0. Hence there exists a real number (eigenvalue) 

kuikywy/ 

{\+\vk) 



A > such that F{\, ko) = 1 and q = ^^'"{^^'y/^ is the eigenfunction. 



We now fix fc = A;o and return to (j5.4p . 

It can be proved (see for example ^) that ||Lg||j;^2 < C||i/g||j;^2. Moreover, for hard spheres, 
there are constants Ci and C2 > such that 

W^^qWl^ = Wikovq- PkoiU{ko){ q\^dC}v^/Jl\\L2 

> Ci\\uq\\L2-C2\\q\\L2. (5.8) 
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This imphes ||Lq||^2 < CIHT^^gll^.a + Ugllia} so the perturbation L is T^-bounded. Since 

T"" =T^ - aL, 

we deduce from Kato's book (P206, [12j) that, for a smaU, there is an eigenvalue A and an 
eigenfunction q{v,C) ^or (j5.4p with positive real part A for ([57 



Theorem 5.2 Let (3 > 1. Then there is a ^-periodic eigenvector {g,g) for —C such that 
ReA > and g{x, v, () = g{—x, —v, C). 

Proof. We fix ko as in Lemma l5.11 The idea is to define, for the family of equations (j5.2p 

ao = sup < a : there is an eigenvalue with positive real part for 

. , 27r . . 1 

with a periodic eigenvector >. 

ko J 



By Lemma [5 .11 for a sufficiently small this set is not empty. We want to show that ao = +oo. 
We prove it by contradiction. Suppose oq < +oo. We claim that, if there is such a finite 



ao > 0, then there is an eigenvalue Aq with an eigenfunction qq with \\qo\\u = y / ^iQol'^d^ = Ij 
such that ReAo = and 

(Ao + ivko)qo - (3kQ\U (ko) | j ^ qo^/JJ'di^ = agLgo- (5.9) 

Proof of the claim: In fact, by (jS.Sp . choose a family of eigenfunctions qa G such that 
\\Qa\\u = 1> RcAq, > 0, as a ^ ao and 

(Aa + ivko)qa - PkoiUiko) | j ^ qay/Jj'd^^ Vy/Ji = aLq^. (5.10) 

Notice that both {Lqa,qa) and (i vkoqa,qa) are bounded by C||(jrQ,||^. As a ^ ao, taking 
inner product with qa for (|5.10p . we deduce that |Aq,| is bounded for a ^ uq. Hence 
linia^ao = Aq (up to subscquences) with ReAo ^ 0- We now prove that Ao is an eigenvalue 
so that ReAo = by the definition of ao and the claim is proven. Clearly, we may assume that 
lim qa = qo weakly in and (j5.9p are valid as a ^ ao . We only need to show that lim q^ = qo 
strongly so that \\qo\\i^ = 1, and qo is an eigenfunction. Denote P/i = {l,f, |vp}y^. Clearly 
Pqa Pqo strongly in L^. It thus is left to show that (I — P)^^ — > (I — P)^o strongly in L^. 
We subtract ( [STOP from (^M) to get 

(Aq, - Xo)qa + (a - ao)Lqa + (Aq + ikov){qa - qo) 
-(3kQ\(j{ko) jy^i^a ~ 9o}\/Atdc| Vy/Ji + aoL{qa - qo) 
= 0. 

We take the inner product with qa — qo and then take the real part. Since ReAoH^'a — Q'olP > 
0, and \ko J v\qa — qol'^d^ is purely imaginary, we obtain: 

"0(^(5" -5), (5" -5)) < (|Aq - Aol + |a - aoDllgalli. • ll^a - golU 

+C\ [ {qa-qo}VJ^dC\-\\ig''-g)\\. 
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Therefore, (I — F){g°' — g} ^ in and II^H^ = 1 and our claim follows. 

Hence, Aq is purely imaginary. Actually we show that it is 0. To do this we take inner 
product with go in (jS.Op to get 



But by integrating ,^//Ix (j5.9p over ^, we obtain the continuity equation 

Ao / qovViJ-dS, + ko'i / qov^/JldS, = 0. 
Therefore, koi J^3 q^v^d^ = Aq /^s QoV^/JIdS, and 



Aolblli-Z^Wfeo) 



qoVJ^d^ 



+ ao{Lg,g) = 0. (5.11) 



Since U{ko) is real and Aq is purely imaginary, taking the real part of (jS.lip we conclude that 
ao{LqQ, %) = 0. Therefore, qo is a linear combination of the collision invariants y^, and 
^-iid 00-^^90 vanishes in (15. 9p . Now ()5.9p reduces to a pure Vlasov equation, and we 
deduce that 

PikoU{ko){Jqo{0^d(}v^ 
Ao + ivjko 

Since / '7o(0\//^<^? / 0) this is compatible with the condition that go is a combination of 
collision invariants if and only if Aq = 0. Thus go(C) = (iU{kQ)yJJl and hence 

pU{ka) = 1 (5.12) 

which is a contradiction to (15. 



Hence we denote (g(f , C)e'^^, A) a pair of eigenfucntion and eigenvalue for the operator 
T"*^. Note that U{—ko) = [/(/co), and {Lq}{—v,(^) = L{q(—v,C)} from the rotation invariance 
of the collision operator. Letting v — > —v in T, we get 

{A - ivko}q{-v, C) + PkoiU i-ko) I j q{-v, O^d^^ = L{q{-v, ()}. 
This implies that q{—v, C,)e~^^^ would satisfy the same ( 15. 3p with the same A. We therefore 

get 

ff(3^,6 = 9(^,C)e''=°" + g(-t;,C)e-''°" 
= {Reg(v, C) + Reg(— C)} cos k^x — {Img(u, C) — lm.q{—v, C)} sin k^x 
+ i{Reg(f , Q) — Reg(— t;, C)} sin k^x + {q{v^ Q + Img(— C,)} cos /cqx. 



is also an (non-zero) eigenfunction since / q{v,C)d£, = 1. Clearly, g{x,v,() = g{—x,—vX)- 

2tt_ 

ko ' 



Hence, for the corresponding period P = such that g = {g{x,^),g{x,^)) is an eigenvector 



for —C. Our lemma thus follows. 



27 



6 Nonlinear Instability. 

In order to establish the non linear instability, we need several lemmas on the properties of 
the fastest linear growing mode. First of all we need to establish the smoothness and long 
time behavior for the growing mode. Recall the definition of the operator £ (jS.ip . We define 

M = {g= [91,92] G L^\gi{x,v,C) = g2{-x,-v,C)} 
and II • Wl'^m will denote the norm on this set. We have the following lemmas: 

Lemma 6.1 Let (3 > 1. Then for sufficiently small, for all 6 > 0, the spectrum of —C 
in {ReA > 5} consists of a finite number of eigenvalues of finite multiplicity. If Ai denotes 
an eigenvalue with maximal real part, and A > max{0, ReAi}, then there exists C\ > such 
that, for any go £ M., 

l|e~*^5o||L2(x) < C'Ae'^*||5'o||L2(X)- 

Proof. This follows easily from the Vidav's Lemma [16]. Notice that we can split 

Cg = {vd^g + Lg} + {PF{^g)v^} 
= Ag + ICg. 

where /C is a compact operator from to L^, while ||e~*'^||j^2(_v()_»L2(_v() < 1- ■ 

Lemma 6.2 Let R = (i?i,i?2) G L'^{Ai) with ||-R||l2(;\^) = 1 be an eigenvector of —C with 
ReA > 0. Then there exists a constant C depending only on A such that 

W'^ccvRWlHm) < C. (6.1) 
supu'(t;)|i?(x,t;)| < C, (6.2) 

x,v 

where w is a polynomial weight as in previous section. 
Proof. We begin with R £ L^. We first claim 

/•oo 

R = - e~^^e-^^ICRdt. (6.3) 
Jo 

Notice that the corresponding growing mode g{t) = e^^R satisfies 

dtg + Ag= -ICg 

so that e^*i? = e-'^^-'^-^e^'R - e'^^-'^^-^ICRdT. Letting s -00, since ||e-*-^||i2^j^2 < 1 
for any t > and ReA > 0, we get 

/t fOO 
-00 JO 

Dividing by e^* we prove our claim because ReA > and the integral converges in L^ . 

From the property of linear Boltzmann equation, clearly 
\\dx{e~^-^g}\\L-2(^M) < WdxgWh'^iM)- Taking v derivative of {dt + vda^ + L)g = yields: 

{dt + vdx]{dyg] + dy{-Lg] = -d^g. 
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Prom [6], / dy{—Lg}dvg > J ^{d^gl"^ — C^\\g\\'^j^2- We thus obtain by taking inner product 
with dyg, 

\\dvg\y = \\dAe-'^9o}\\L^ <C{t + l){\\go\\L2 + \\V,,,go\\L2}. (6.4) 

Since ICR £ C°° and ||52:{/Ci?}||j;^2 + ||(?t,{/Ci?}||j;^2 < C||i?||/^2, we can take dx and dy derivatives 
in (|6.3|) to get 

/•oo 

\\dxR\\L2 + \\d.R\\L^ < / e-^'^'{\\dx{e-'-^ICR}\\L2 + \\d,{e-'-^!CR}\\L2}dt 

Jo 

< C / e-^^^*(t + l)\\dx{lCR}\\L2 + \\d,{lCR}\\L2dt 

Jo 

POO 

< C\\R\\l2 e-^^^\t + l)dt<C\\ R\\l2. 

Jo 



We therefore deduce (j6.ip . To show (16. 2p . we denote S = wR. We then have 

\S = {vd^S + wL{-)} + {pF(^^)vwJJl} = A^S + ICyjS. 

w w 

Applying the same proof in Section 3 for the stabiUty for the pure hnear Boltzmann operator, 
we can estabhsh: 

1 1 e"*-^"' 50 1 1 L- < <^{||5oIIl2 + Ibolli-}- 
We can similarly obtain S = — e~'^^e~^-^^ ICy^Sdt, so that from ReA > 0, 

/•OD 

Jo 

POO 

< C e-^'^'{\\IC^S\\L^ + \\IC^S\\L2}dt 

Jo 

< C. 



Lemma 6.3 Let R be an eigenvector of —C with its eigenvalue A with ReA > 0. If \ is not 
real, then there is a constant C > such that for all t > 0, 

\\e-^HmR\\L2 > Ce^^^*||Imi?||2 > 0. 

Proof. We prove by contradiction. Notice that, since 
lmF{yJJIR) = F(y^Imii), one can immediately check that 

e'^HmR = lm{e^^*R} = e^''^*(sin[ImAt]Rei? + cos[ImAt]Imi?). 

If the Lemma were false, by passing through a convergent subsequence of 

sin[ImAt„] and cos[ImAt„], with n — > oo we would have almi? + 6Rei? = 0, with + 6^ = 1. 

Therefore either ImR or Rei? would be a real eigenvector and A would be real, a contradiction. 

■ 

Lemma 6.4 Let R be as in the preceeding lemma. There exists Sq > 0, such that for < 6 < 
6o, there exists a (compactly supported) approximate eigenfunction Rs such that 

S\Rs{x,v)\.y]2 < fi, 

WR-RsWl^ < V6, 
\\dM\L^ + \\dM\L2 < C{\\dxR\\L^ + \\d,R\\L2}. 
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Proof. In fact, we choose x(^) be a smooth cutoff function x(^) = 1 for \v\ < N and xi'^) = 
for \v\ > N + 1. By Lemma 16.21 we have 

\x{v)R{x, v)\^< \Rix, v)\^= \wS\^ < {Cwfi-'/^}ti. (6.5) 

Define by the equation 6 = ^q^^^I^^I^ and define 

Rs = x{v)R{x,v). 

Clearly the third estimate in the Lemma is valid. From ()6.5p and the definition of N and 
S, the first inequality in the Lemma is also valid. Since is a polynomial, we have y/6 = 

\l ^Cw{N^ ^ — V l^i-^) ■^liGii ^ is large. We then conclude the lemma by 

\\R - Rs\\l2 = ||iil|,|>^||L2 <C [ ^^^dv = C^i^iN) < Vs. 

J\v\>N W{v) 

m 

We now establish the crucial bootstrap Lemma which shows that growth leads to the 
same growth rate for L°°. 

Lemma 6.5 Let g = (51,32) be a solution to the nonlinear problem around Mhom-' 

(dt + vd^Qi + PF{^gi+i)^v - Ug (6.6) 
= -F{VJ^9i+i)dv9i + F{^gi+i)giV + T{gi,gi) + r{gi,gi+i). 



Assume that ReA > and 

II _/_i\ii ^ ^_ReAtii _/n\i|^^ 



\g{t)U.<Ce''''% 



fort G [0,r]. There exists Eq > such that i/ supo<t<j'{||it;(?(t)||Loo + ||5(i)||L2} < eq, then 
there is a constant C such that 

\\dx9{t)\\L2 + \\dvg{t)\\L2 + \\wg{t)\\L^ 

< Ce^'^'mgrniL^ + 115,5(0)11^2 + ||/i(0)|Uoo}. (6.7) 

Proof. We take x and v derivatives for ()6.6p . Since F{M.i) = 0, from (j4.16p . we have {h = wg) 
^^Wd.gWh < C\\h\\L^{\\d,g\\l, + \\d,g\\l,} + C\\g\\l,, (6.8) 
J^WdvgWh < C\\d^9\\h+C\\g\\l,. 

Applying the Gronwall's inequality to (j6.8p . by ||/i||loo < eg < ReA, we obtain 

t 



Id^ml, < \\dMO)\\l, + Ceo I e^^«(*-^)||a,5(s)|li2ds 

I /^„ReAt II „/nM|^^ 







+ Ce^-^^\\9m\l„ (6-9) 
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+ \\d.9m\h}- (6-10) 

Letting Y{t) = £ e^^o^*"^) ||5^,g(s)|||2(is so that Y' = \\d^g{s)\\l2. Multiplying with Y' > 0, 
we have 

Gronwall's inequality and integration by parts implies 

{Y'}\t) < \\dMO)\\l.+C{\\gmh + \\dMO)\\hV Te^^^^'-^^e^^^^^d. 

JO 

< CmgmWh + |b(0)||i. + 115.5(0) Hi. Pe^^^^*. 

Taking the square root and plugging into (|6.9|) , we thus conclude the desired estimate for 
\\d,g{s)\\l, + \\dMt)\\l, in (I62D. 

We now consider the L°° estimate. By repeating the same (but simpler) argument in 
Lemma [4. II for the full nonlinear system (with F{Mi) = 0) , we can show that if {||/i(t)||L°o + 
115(011^2} is sufficiently small, 

\\h{t)\\L^ < C{\\hm\L^ + sup \\g{s)\\L2} < Ce^'-^'mmL^ + lb(0)||L2}. 

0<s<t 

■ 

Proof of Theorem 11.41 We choose R to be the eigenfunction whose eigenvalue has the 
largest positive real part. If A is not real, then ||Imi?||2,2 = r > 0. We choose the approximate 
eigenfunction ImiZ'' to the imaginary part of R by Lemma 16.41 In case A is real, we simply 
do not take the imaginary parts. 

We choose a family of solutions /''(O, j;,^) = /x + 51m.R^ ^JJi > or g^{0,x,^) = 6lmR^ . 
Note that the positivity follows from the first statement in Lemma [6.41 for 6 sufficiently small. 
Clearly, from Lemma 16.41 

11/(0) - 6lmR\\L2 = 6\\R - Imii-^H^. < 6^+^ < y , 
for 6 sufficiently small. Hence, from Lemma 16.41 

||/(0)lbi + 1 1 /i"^ (0)1 1 LOO = 6\\lmR^\\Hi + 6\\lmR^\\Loc < C6r. 



Now from the nonlinear Vlasov-Boltzmann system ()6.6p . we have 

g^{t)=6e-^HmR^ (6.11) 

^-C{t~-T) (-F{y/J^92)dvgi + F{^g2)vgi+T{gi,gi) + T{gi,g2)\ 
Jo \-F{^gi)dvg2 + F{^gi)vg2 + T{g2,g2) +T{g2,gi)J 



+ 

We choose A such that 
Let 



ReA < A < (1 + ^)ReA. (6.12) 

r'5^ = ^4^lin. ' 



A -ReA' 2CaVS 
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By (I6l^ and = ^In f . Since 2(A-ReA) < ReA, for small 6, we have < T^^ . Clearly, 

g(A-ReA)t < g(A-ReA)T*^ ^ for < t < T^^: 

2Ca V o 

Ca^/^c'^* < ^e^^^V. (6.13) 

Let 

T* = sup{s : ||V,/(t)||i2 + ||V„/(t)||i2 + Wh^mL^ < so} (6.14) 

s 

T** = sup{s : \\g^{t) - 6e-^^R^\\L2 < j5e^^^^r, for all < t < s}. (6.15) 

For < t < min{T^r**}, we have from (l6l^ . 

||/(t)||L2 < <5||e-^*Imflii2 + ^5e^^^V (6.16) 
= 6\\e-^HmR\\L2 + 6\\e-^'{R - R^}\\l2 + ^<5e^"^V 



< (l + ^)ei^^^*||/(0)||^.. 



4 

We now claim that < min{r*,r**}. 

In fact, if T* < minfT-^, T**}, then by (fOell . (I6l^ . (f636|) . and the Lemma[63] (bootstrap 
lemma), we obtain for < i < T*: 

l|V.,./(t)||L2 + \\h'{t)\\L^ < C7eR<=^*{l|V.,./(0)||i2 + \\h\0)\\L^}. 

In particular, by p.26p . 

\\y.,v9'iT*)U2 + ||/i'(r*)||ioo < Ce^^^^*{||v,,,/(o)|U2 + \\h'm\L^} 

< Ce^eAT^5 = C0. 

This is a contradiction to the definition of T* if is chosen « eq. 

On the other hand, if T** < min{r'^,r*}, then by (l6l^ . (I6T^ . and LemmaESI 

||/(r**) - 5e-^^"lmi?^||i2 < C {\\V,,,g\t)\\L2 + \\h\t)\\Loo}^dt 

Jo 

< Ce2R-^^"{||V,,,/(0)||i2 + Wh'miLoo}^ 

This is a contradiction to T** in (j6.15p when 6 is small. 

Now that < min{r*, T**}, we can evaluate t = T"^ in (f02]l to get 

||/(T^) - 5e-^^'lmi2''||i2 < Ce\ 
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But 

\\6e-^^^'"lmR^\\L2 > \\5e-^'^" lmR\\L2 - \\Se-^^'lm{R - Rs}\\l2 

4 4 

Therefore, for 9 sufficiently small, 

\\g\T')h2>^re-ce^>^re. 
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